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LIS 1 Ob BUJLNCll’LL SYMBOLS 


n = number of buses in the s>stem 

& V = voltage angle, magnitude at bus 1 

e^ , f i = real and imaginary components of the voltage 

Y = (j + jB = ’ii th element of bus admittance matrix 

ij ij ij 

1* + jQ^ = net real and reactive power injected at bus l 

p + jq = real and reactive line power flow from ith bus to 

i J ij 

jth bus 

lhe other symbols in the text are defined as when the> are 
introduced 
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Ihe energy control center is the nodal agency which controls 
the power system network lhe load flow solution program is run to 
determine the power flow in the network lhe input to this 
algorithm is the measured data telemetered from field and may get 
corrupted due to various reasons Hence, the state estimator 
program is used to estimate the state of the network Ihus, the 
load flow and state estimation are two subproblems of power system 
network 

In this work, the performance of various load flow and state 
estimator algorithms are critically compared in well and 
ill-conditioned cases lhe standard fast decoupled load flow in 
polar coordinates using Stott’s as well as Amerongen s assumptions 
are compared with the ^DLl* methods in rectangular and two new 
algorithms m hybrid coordinates 

lhe literature survey shows that orthogonal isat ion by Givens 
rotation has been applied to the state estimation using normal 
equations tast decoupled state estimator by normal equations 
method using the same assumptions of Stott as well as Amerongen is 
compared with the normal equations using orthogona l isat ion 

A comparative study of all the bULl* and the various state 
estimator methods has been done in this thesis with reference to 
computation time and number of iterations for both normal and ill 


conditioned cases 
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IN i ltODUC 1 ION 

1 1 GENERAL 

Ihe electrical power system is undergoing continuous 
expansion m size and complexity Power system networks are 
being provided with computer aided hnergy Management system 
for its effective monitoring and control it is essential that 
the operating state of the system be continuously monitored, so 
that this knowledge can be be used for operation and control of 
the system Energy control centers receive measured data 
from field through remote terminal units and communication link 
Determination of system state in absence of any measurement 
error can be done by load flow or power flow studies whereas it 
can be performed by state estimator m presence of random 
error m the measurement 

Since the operator’s decision for planning, operation and 
control are based on the knowledge about system operating 
conditions and many other advanced functions are performed 
based on the knowledge of current system state, the two 
subproblems viz, the load flow and state estimation are of 
vital importance lhe present thesis addresses to the two 
problems subproblems and explores few fast methods suitable 
for their on line applications 

Load flow study is a steady state solution of a power system 
problem which provides information about network power flow and 
voltage profile Load flow studies on a digital computer 
essentially involve the solution of non-linear algebraic power 
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flow equations It is the most frequently carried out study for 
taking various on-line as well as off-line decisions 

State estimation is a software, run on-line at energy 
control centers and is a process of assigning values to the 
un-known system state variables using imperfect measurements 
the state estimator requires the measurements to be redundant 
and the process of estimating the system states is usually based 
on a statistical criteria that estimates the true value of the 
state variables to minimizes or maximize the selected 
criteria A commonly used and familiar criteria is that of 
minimizing the sum of the squares of the differences between 
the estimated and measured values of a function 

In the power system, the state variables are the voltage 
magnitudes and phase angles at the system buses the inputs to 
the estimator are selected imperfect power system measurements of 
voltage magnitudes, real power, reactive power or ampere-flow 
quantities lhe state estimation can smooth out small random 
errors in measurements, detect and identity gross measurement 
errors and fill-in meter readings that have failed due to 
communication breakdown 

lhe desirable features of these two subproblems are 
1) High speed for real time application 
u) Low and compact storage 

111) Reliability for well behaved and ill-conditioned situations 


m the networks 
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1 2 STATE OF THE ART AND MOTIVATION 

in the past four decades, a number of research works have 
been reported on load flow and state estimation problems it is 
not possible to review all the work reported m these two areas 
Hence only a representative survey of few important works and 
those relevent to the present thesis are presented below 
121 LOAD FLOW 

lhe history of load flow analysis dates back to the fifties 
when method, now known as Lauss-Seidel [1] was used for solving 
the load flow equations Although this method was very simple and 
requires very little memory storage, it requires large 
computational time and presents convergence problems for large as 
well as ill-conditioned systems 

lhe Newton-ltaphson method represented a breakthrough [31 in 
load flow situation lhe Newton Uaphson load flow (.NULl) method 
involves the linearization of the set of load flow equations at 
any iteration point using only the first order term m laylor's 
series expansion Several versions of NH method have been 
developed and the popular ones are based on power mismatches 
lhe power mismatches versions can be formulated using either 
the polar [3] or rectangular co-ordinates [21 for voltages lhe 
power mismatch version in polar co-ordinates has advantage over 
rectangular version in that it requires lesser memory space 
as the number of variables are reduced 

lhe power mismatch version m polar co-ordinates developed 
by linnej and Hart [31 utilizes the sparsity exploitation of 
optimal ordering of buses and has been widely used by the 
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utilities lhis method usually converges to a fairly accurate 
solution in two to five iterations for most of the systems 
independent of their size and type though NULF method has 
proved to be highly reliable, the storage and time 
requirements were considerably high 

Decoupled versions of NU method t5,b] were proposed in view 
of reducing memory and computational time in all these methods, 
the weak coupling between the active power - voltage magnitude 
and reactive power - voltage phase angles were exploited Amongst 
various decoupled versions, fast decoupled load flow ODU 1 ) 
method developed by Stott and Alsac [b] is the most popular with 
utilities Ihe method is m polar co-ordinates and utilizes 
certain justifiable network assumptions apart from the l'-V and 
Q-<5 decoupling to obtain constant Jacobian submatrices lhe 
final model of tBLi* is developed employing certain additional 
assumptions m computing elements of the submatrices [B*] and 
[B’’J lhe frDLl* model developed is quite fast and reliable for 
a wide spectrum of systems 

Over the period from 1975 till date, a number of 
modifications have been proposed to Stott *s fast decoupled load 
flow to improve its speed, convergence characteristics and 
memory requirements in addition to this, a number of new 
but similar decoupled versions have been proposed which easily 
lend them-selves to comparison with Stott’s IDLl* 

Nanda, et al [111 have in their paper brought out clearly 
the relative importance of some of the assumptions made by Stott 
and Alsac in their tDLfr lheir paper states that the omission of 



shunt reactances and off-nominal in phase transformer taps m 
[B ] matrix has little effect on convergence of the 1 UL1* and the 
omission of series resistances, however, plays a vital role m 
reducing the number of iterations for certain systems 

Van Amerongen [18] has studied the effect of neglecting the 
line resistances and has come up with four schemes wherein the 
line resistances are neglected or included in either or both the 
[B’ ] and [B’’l matrices lhe test results show that for higher 
H/X. ratios the scheme where the line resistances are only 
neglected m [B '] matrix perform better than as suggested by 
Stott and Alsac lie also advised the strict Cld-lV) iteration 
scheme to avoid cycling behavior lhe paper also suggests tbu~+ 
cycling can be avoided even with the standard iteration schc’ne 
suggested by Stott and Alsac, by using different tolerances for 
the active and reactive subproblems whii.-* carrying out 
iterat ions 

Hubbi [21] has theoretically s. own the differences between 
the Stott’s and Amerongen’s H)L1* Keyhani [15] has come up with 
four modified versions ■‘f the MIL! 1 in order to reduce the memory 
requirement for very large systems lhe algorithm is claimed to 
require „ less memory as compared to the HJLfr However, it has 
a slightly slower rate of convergence than the IDLl 

An improved version of the frDLl- by Mark knns [8] presents a 
new algorithm based on a slight modification of the MlLk lhis 
eliminates the calculation of trigonometric functions altogether 
and yet converges to an exact solution lhe convergence property 
is same as for the ^DL* and the savings in multiplications is 
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about l/4th for each iteration thus computation time is reduced 

Behnam and (lUliani tlbl have described a load flow method at 
least 50% faster than btott and A1 sac’s M)Lt' lhe model used by 
the authors combines the nodal iterative model (Gauss— Seidal 
technique) and the tDLt* by btott and Alsac An important and 
useful feature of this algorithm is that this model does not use 
[B 1 matrix and therefore, m contingency studies, the [B’*l 
matrix needs to be refactor lsed for accommodating the 
contingencies, this method gives a faster solution than the 
standard (<1)1,1* 

lhe fast decoupled load flow of Stott and Alsac is based on 
the assumption that the U/X ratios are small for all the 
branches lhis condition constitutes a limitation m finding load 
flow solution for those which have branches having relatively 
high resistances or where the overall 11/X ratio is not small 
lhis limitation has elicited many research works m this area 

Van Amerongen [18] has come up with one solution as stated 
earlier Dy Liacco et al [7] have suggested the series 
compensation method to overcome the problem, but could not 
compete over the speed {similarly, tlonticelli et al [10] have 
come up with shunt compensation method, which is more reliable 
than the series compensation method llajicic and Bose [17] have 
presented a new modification to the tDLt , which they found, 
consistently provided better convergence Iheir modification is 
based on tests conduct-ed by the authors on systems with high lt/X 


ratios However, their method converges to the solution only for 
U/X ratio below three Wang et al [22] too have developed 



similar technique on further modifying the algorithm proposed by 
Anjan Bose et ai [15] 

from the literature survej it is found that the fDLf method 
based on Stott and Alsac s assumptions [b] are still most popular 
compared to the other methods and is being used for various 
on-line and off-line studies lhe General purpose fDLf method 
suggested by Amerongen [18] offers advantages m solving 
systems with high ll/X ratio of lines and is also being adopted by 
utilities few models suggested in literature do not consider 
Q-limits of generators 

Although much research efforts have been expended m 
developing fast decoupled load flow models in polar 
coordinates! the development of its counterpart in rectangular 
coordinates has not attracted much attention 

lhe mam advantage of the rectangular version over the polar 
version is that the former does not involve trigonometric 
functions and thus leads to a considerable saving in time spent 
on the computation of these If decoupling can be realized 
effectively, then the rectangular fDLf versions would be as 
competitive as the polar fDLf versions 

llaju et al [14] have developed the first order decoupled 
algorithm whose model is based on all the assumptions made by 
Stott lheir algorithm is in rectangular co-ordinates but the 
scope of their algorithm is limited to well-behaved systems 


Bab ic 

[13] 

too has 

come up with decoupled 

load flow 

m 

rectangular 

co-ordinates , 

which however is 

not 

re 1 table 

for 

1 1 1-condit toned 

systems 

L Srivastava et 

al 

[20] 

have 



formulated three 


different fast decoupled algorithms m 
rectangular co-ordinates, one of which exhibits better 
convergence and is superior to the Stott's IDLI in certain 
ill-conditioned cases in terms of number of iterations but 
takes more CPU time Also the methods did not consider Cj- limit 
of generators 

Since the power flow equations in rectangular co-ordinates 
involve only quadratic terms, the complete laylor’s series 
expansion will be limited to only second order terms, unlike the 
polar co-ordinates hence the second order load flow m 
rectangular co-ordinates, m effect should be more accurate than 
the polar or first order load flow m rectangular 
co-ordinates Some if the decoupled versions of second order 
load flow in rectangular coordinates were suggested by Nanda 
et al [121 and Cory et al [19], but is not found to be 
computationally superior to the IDLI in terms of Cl’U time 

lhe main advantage of the rectangular version over the polar 
version is that the former does not involve any trigonometric 
functions and thus leads to a considerable saving in time spent 
on the computation of these If decoupling can be realized 
effectively, then the rectangular IDLI versions would be as 
competitive as the polar IDLI versions 

Hence the motivation was to develop two IDLI models m 
rectangular coordinates using Stott s as well Amerongen’s 
assumptions and considering generator Q-l units In addition, two 
hybrid algorithms where in the two decoupled equations are in 
two different coordinates have been developed lhe performance 



of all methods have been compared with the s-tandard tDLl- [b,18] 

in polar coordinates 

12 2 STATIC STATE ESTIMATION 

Ihe static state estimation (SSL) considers the time 
invariant measurements and static model of the power system One 
of the earliest methods of state estimation algorithm 
developed by Schweppe et at [23] is the Weighted Least Squares 
(WLS) method lhe WLS algorithm works well and converges fast 
to a solution point with the systems whose information matrix 
is well conditioned When the information matrix is 
numerically ill-conditioned, the WLS algorithm does not 

perform satisfactorily and manifests this problem m the poor of 
slow convergence 

lhe fast decoupled static state estimator developed by 


Horisberger 

et al 

[24] 

utilizes the decoupled l’-d 

and 

Q-V 

equations 

and 

solves 

the 

active and reactive 

equat ions 

iteratively 

using 

constant 

simplified sub-matrices 

of 

the 

inf ormat ion 

matrix 

lhe 

algorithm yields the exact solution 

and 


it requires less computer storage However it could not detect 
and identify bad data 

Garcia [25] introduced fast decoupled state estimation using 
decoupled detection and bad data identification Bad data is 
eliminated by pseudo measurement generation and avoids gam 
matrix retnangulation 

llao et al [2b] described an improved decoupled static state 
estimator based on minimization of WLS residuals No 


approximation is made 


in the mismatch 


functions 


the 


f mal 
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solution is as accurate as provided by using complete 
information matrix lhe estimator requires considerably less 
storage and solution time but more iterations Ihey have shown 
that the rate of convergence of the estimator is significantly 
improved b> adaptively adjusting the step-size during 
iterat ions 

Monticell 1 et al [33] developed fast decoupled state 
estimator in a new frame of reference Decoupling was not 
seen as zeroing coupling submatrices in the problem Jacobian, 
but as a two step procedure to solve full Newton equations, 
without major approximations in the primal algorithm instead 
of neglecting resistances from [B’ ] matrices, resistance is 
neglected from [B*’] matrix Whereas dual algorithm is similar 
to standard fast decoupled state estimator 

bnnxvasan et al [31] developed three new and improved 
methods for static state estimation problem, out of which two are 
based on first order model and the third one is based on second 
order model Among the three methods the modified iast Decoupled 
State hstimator is the fastest algorithm though the other two 
methods are not as fast as the modified t*DS estimator, they are 
more reliable than the fcDS estimator, l e , well suited for 
ill-conditioned system 

Irving et al [25] introduced power system state estimation 
using linear programming Linear programming has the combined 
advantages of noise filtering and bad data elimination and was 
implemented using the Simplex method Ll-Keib et al [37] came up 


with 


a 


formulation which improved the 


solution 


t me 
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significantly Ihe disadvantages of the LP estimator are that, 
it does not exploit the problem structure and it fails to reject 
the bad data due to presence of leveraged points 

Kotigua et al [301 used a new method for obtaining an 
estimate of the state of a power system using Weighted Least 
Absolute Value tWLAV) WLAV estimator simultaneously detects and 
rejects bad data while obtaining an accurate estimate of the 
state Celik et al [35] developed a robust state estimator using 
transformations, which remains insensitive to bad measurements 
even when associated with leveraged points these 
transformations represent a change of co-ordinates m the 
state-space WLAV estimators lose their robustness m the 
presence of leveraged points lhe transformation of Jacobian 
matrix [H] is to eliminate the leveraged points associated 
with the original set of measurements WLAV estimator’s draw 
back is its poor computational efficiency for large systems 

boliman et al [3b] described a new technique for constrained 
power system state estimation by considering equality 
constraints Ihis is based on WLAV state estimation procedure 

llao et al [38] described the Levenberg-Marquardt algorithm 
to determine the static state estimation for i ll-condit loned 
power systems lhe diagonal elements of the information matrix of 
WLS estimator are suitably increased to accelerate the 
convergence during the iterative solution when the linear 
model is found to be too approximate 

Van Lutsen et al [39] used a hierarchical concept to 
solve the static state estimation problem for large scale 
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composite power systems In this algorithm solution is obtained 
by performing a two level calculation In the lower level, a 
conventional state estimator is carried out simultaneously 
for all sub-systems lhe co-ordination of these local 
estimations is realized in the upper level, thus it receives 
and treats only a small number of variables 

Labudda et al [38] used fuzzy multi-objective approach to 
power system state estimation lhe fuzzy LI* estimator offers a 
number of advantages over the standard Ll’/WLAV approach, the most 
noticeable of which is superior robustness in the presence of 
leveraged bad data 

LarsHolten et ai [32] compared different methods for state 
estimation such as Normal equations method, Orthogonal 
transformation method, hybrid method, Normal equations with 
constraints and Hachtel’s augmented matrix method 

lhe normal equations method of solving the static state 
estimation does not permit the calculation of the power system 
state variables when new measurements are obtained after the base 
case has been solved ihe equations technique requires a complete 
rerun of the solution algorithm with regard to the new 
measurements 

However, the orthogonal transformation method permits the 
recalculation of base case state variables with respect to the 
new measurements Narasimhan Vempati et al [34] have developed 
the transformation using the Givens rotations lhey have 
suggested three enhancements, to improve the speed of 
orthogonal isat ion for power system state estimation, and make 



11 comparable with that of the normal equations method 

Iron the literature survey, it appears that the PDSE is 
being popularly used by utilities Amongst the recently developed 
methods, the orthogonal transformation method [34] shows promise 
in achieving the solution much faster than other methods 
However, the model suggested in the paper was based on Nil 
method Hence the motivation was to extend the orthogonal 
method using frast Decoupled equations using both Amerongen’s 
and Stott's assumptions and compare its performance with 
ibSJb based on normal equations 
1 3 THESIS ORGANIZATION 

lhe work reported in this thesis has been organized into 
four chapters 

lhe present chapter - 1 introduces the problem of load flow 
studies, state estimation and sets the motivation behind the 
present work 

Chapter - 2 is devoted to the comparison of different last 
Decoupled Load How models two algorithms for 1DL1 m 
rectangular coordinates and the two using hybrid model have 
been developed m this chapter Potential of all these 

methods have been tested on IEEE 14 Bus and 30 Bus test systems 
for both well-behaved and ill-conditioned cases 

Chapter - 3 deals with the last Decoupled Power bystem State 
Estimator two different Models of static state estimator has 
been proposed using orthogonal transformations based on Givens 
rotations and its performance has been compared with the PDSE 
using normal equations method lhe tests have been carried on 
IEEE 14 Bus test system 
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frinall>, the Chapter - 4 concludes the specific findings in 
this thesis along with suggestions for future scope of work 
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COMl’AlliSON Ob DiibLHLNl 1 JDLt MklHODS 

Z 1 IN 1 HODUC i ION 

lhe objective of load flow stud* is to determine the voltage 
magnitude and angle at each bus, real and reactive power flow in 
each line and line losses of the power sjsten network for 
specified bus loading conditions lhe Newton llaphson's method has 
proved to be a powerful method for solving the load flow equations 
but is not used for on-line applications as it requires more 
memory and computational time 

bast Decoupled Load blow IbDLb) method suggested by Stott and 
At sac [b] overcomes these drawbacks of Newton Haphson method and 
is probably the most popular method being used by utilities 

this method utilises certain physical properties of power 
system networks such as 1 ’— & and Q— V decoupling, considers 
assumptions such as lines have small H/X ratio, voltage angle 
difference at two buses being small In addition the model 
utilises certain additional assumptions to make the Jacobian 
matrices [B'J and [B ] constant for which no theoritical proof is 
given by Stott and Alsac [b] Several other models considering 
these additional assumptions m different combination were 
suggested and the one become popular is a General purpose bDLb due 
to Amerongen et al [18], which was subsequently proved by 
Monticelli et ai l J 

lhese models were developed in polar coordinates and use 
sensitivit* method to handle Q limits of generators However, 
their counterparts in the rectangular coordinates considering Q 


limits have not been tried at hence, in this chapter the IDL* 
models based on Stott's and Amerongen's assumptions have been 
developed m rectangular coordinates and tested for Ihfcfc. 14 and 30 
bus systems lwo hybrid models IDLl models using both polar and 
rectangular version of l*-«5 and Q-V equations have also been tried 
and their performance compared with the above IDLI methods in 
polar and rectangular coordinates 
2 A IDUb Mil 1 HODS IN 1‘OLAH COOltDINAlhS 

In a practical power system, the following assumptions hold 

good 

i) (6 - 6 ) is small, thus C0SC6 - 6 ) 3? 1 

t J V J 

ll) C S1NC6 - 6 ) « B 

IJ t J M 

ill) Q «B*V 2 

L LL L 

Consider a power network system to have a total n+1 buses 
with the (n+l)th bus as the slack bus, Cl-m) being l’-V type and m 
I»Q type buses lhe following equations in polar coordinates are 
solved to get the unknown variables l e the voltage magnitudes and 
phase angles 


‘ Al> 



r 1 


V 


M 


C2 1) 


' AQ ‘ 





V 

— * 

[“ ■] 

H 

12 2) 


In the above eqns [B* ] and [B* ’ 1 matrices are constant, real 
and sparse matrices and need to be triangulated/ inverted only once 
at the beginning of the iteration 

Stott and Alsac [b] took the following additional assumptions 


to arrive at the above mentioned equations 
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l) Omitting from [B ] the representation of those network 
elements that affects MV A11 flow 1 e, shunt reactances and off 
nominal in-phase transformer taps 

n) Omitting from [B ’1, the angle shifting effects of phase 
shifters 

m) Neglecting series resistances in calculating the elements of 
IB’ 1 

lhe above equations form the fast decoupled load flow model 
[&] on polar coordinates hach iteration cycle comprises of one 
solution for [A6] to update [6J and then one solution for lAV] to 
update [V] and is termed as (16, IV) scheme lhe iterative 
procedure is repeated till the real and reactive power mismatches 
become less than a pre-ass igned tolerance value 

Stott and Alsac’s model of kDLk however could not perform 
well m case the power system network was ill-conditioned, l e , 
with high H/X. ratio Van Amerongen [181 has formulated a general 
purpose kDLk algorithm in which he has neglected the series line 
resistance from the IB' ] matrix and not from [B*l matrix 
Amerongen’ s model of kULk was able to solve the load flow problem 
even m ill-conditioned cases lo compare the various kDLk models 
critically for their performance m well and ill conditioned cases 
the Stott's model and Amerongen’ s model are termed as kDLk-i and 
kDLk-11 respectively 

2 3 kDLk ML 1 HODS IN HLCIANOULAU COOUD1NA1LS 

kast decoupled load flow method [b,18j in polar coordinates, 
kDLl‘, discussed earlier is found to be very popular method because 
of their less memory and Cl’U time requirement lhe formulation of 
a load flow method m rectangular coordinates offers faster 
solution as compared to its formulation m polar coordinates, 
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mainly because of the reasons that it does not involve additional 
computational efforts in solving trigonometric functions With 
this in view, development of fast decoupled load flow methods in 
rectangular coordinates OULU) have been tried out 

bBLl* model in rectangular coordinates can be developed using 
the similar approach as in polar coordinates Ihe Voltage at any 
bus m rectangular coordinates can be considered as 

V = e + jf <2 J) 

1 v l. 

Assuming that e = 1 0 p u and Q « Be 2 the load flow 

V V VL V 

jacobians are found to be nothing but that derived from negated 
susceptance part of the Y matrix taking additional assumptions 
as done in the standard tDUb m polar coordinates, the t-DLl- model 
m rectangular coordinates are 


M -HM 


(2 4) 


M - [“'■] H 


(2 5) 


lo determine the voltage magnitude of the 1*-V type buses, the 
fast decoupled load flow method m rectangular equations has m 
addition to the eqns 12 4) and 12 5), one more equation which is 
e 2 = tV Bp > 2 - f 2 (2 b) 

lo consider the Q limits of the 1*-V buses, two new models m 
rectangular coordinates using the sensitivity method has been 
proposed in this thesis Ihe two model of bDLb in rectangular 
coordinates taking the Stott’s and Amerongen s assumptions are 
termed as bDLb-Ili and tDLb-lV respectively 


2 4 1‘UOl‘OShl) bASi DhCOUl'LhD MklllODS 


As has been stated earlier, the fast decoupled method in 



lj 


polar coordinates is less faster than the formulation in 
rectangular coordinates because of the additional computational 
efforts in solving the trigonometric functions However, the 
number of iterations needed for convergence m rectangular 
coordinates is found to be greater than that in the polar 
coordinates Hence, it was with a view to combine the advantages 
of both the coordinate method that two new formulations in hybrid 
form have been proposed and has been critically examined for their 
performance in well and ill-conditioned systems 
2 4 1 PROPOSED HYBRID 1-DLE ME 1 HOD 1 


In this proposed method, termed EDLt-V 
mismatch equation is in polar coordinates and the 
mismatch equation is in rectangular coordinates 
are as follow 



M - [»"] H 


the real power 
reactive power 
lhe two equations 


(2 6 ) 

(2 7) 


A flat voltage (1 0 + j 0 ) for P-Q buses and V° poe for slack 
and P-V buses is initially started with Since the eqn (2 b) is m 
polar coordinates, the real power is calculated using bus voltages 
m polar coordinates and the real power mismatch is determined 
After the first half iteration is completed, the incremental bus 
phase angle CA5) value is used to update the complex bus voltage 

lhe complex bus voltage m polar coordinates is converted to 
rectangular coordinates and the reactive power is calculated 
Subsequently, the reactive mismatch is determined and the second 
half iteration is carried out to obtain (Ae) this is continued 
till convergence Major computational steps for the proposed model 



n 

Li 


u 


of tDLt is summarised as fol towing 

1) Head system data including line data, loads at I’-Q buses, 
generator data etc 

2) torm ¥ hos; matrix, [B ] and CB' * ] matrices and I’-V buses 
sensitivity values 

3) initialize the bus voltage to Cl 0+j0 )p u for all l’-Q 

buses, for slack and l’~V buses 

4) Set iteration count K = 0 

5) Calculate [Al>] for all buses except slack bus with bus 
voltage m polar coordinates 

b) lest for convergence of bus real powers (max Al’ < e) If 
converged, go to step 10, else go to step 7 

7) Advance iteration count by 05 ( K = K + 1 ) 

8) Solve for A«5 using eqn <2 b) 

9) Update <5 by A 6 for l’-Q and l’-V buses 

10) Convert the complex bus voltage in polar coordinates to 
rectangular coordinates 

11) Calculate AQ for all buses except slack bus and check for 
convergence If converged goto step lb else go to step 13 

12) tor 1*-V buses, check for violation of Q limits, if so 

determine Ae by using the sensitivity eqn 
AQ 

Ae = * S where S is the sensitivity of the ith i’-V 

L t 

e 

bus 

13) Advance the iteration count by 05 (K = K + 1) and compute Ae 
using equation <2 7) 

14) Update e by Ae for all l’-Q buses 

15) Convert voltages in polar coordinates and go to step 5 

lb) Check for convergence of eqn (2b) If not converged go to 


step 5 
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in this proposed method 1 termed EDLE-V1) the real power 
mismatch equation is m rectangular coordinates and the reactive 
power mismatch equation is in polar coordinates Ihe two mismatch 
equations are 


M ■ [ “'] [ Af ] 


12 9 ) 


AQ ' 

V 


[• ] M 


12 10 ) 


A flat voltage <1 0 +J0) for P-Q buses and V° poe for slack 
and 1*-V buses is initially started with Since the eqn 12 9) is 
m rectangular coordinates, the real power is calculated using bus 
voltages m rectangular coordinates and real power mismatch is 
determined lhe bus voltage is updated by Af after the first 
iteration 

lhe bus voltage in rectangular coordinates is converted to 
polar coordinates and reactive power mismatch is calculated lhe 
second half of the iteration determines the incremental voltage 
magnitude AV 

lhe computational steps for the model iULi-Vl are similar as 
discussed in the ear liar proposed method 

) 

4 


2 b HANDLING OE PV BUSES 



*m. No 


real 

power capacity 

bus , 

and 

the other 

in 

the 

load flow 

for 

the P 

-V buses 

generat ion 

and load 


lhe reactive power limits of the generator are known in the 



process of calculation of the load flow solution 


if 


the 



calculated reactive power violates the specified limits for any 
l'-V bus, then two techniques are employed to overcome this 
problem 

i) Bus switching method 
11) bensitivity method 

Ihe Bus switching method is the technique wherein m each 
iteration the calculated reactive power for the P-V bus which 
violates the limit (.either the maximum or minimum value) is fixed 
to its corresponding limiting value and the l'-V bus is treated as 
a P-Q bus during that iteration and the change m the voltage 
magnitude is determined 

With the bus switching scheme, if during an iteration a P-V 
bus violates the Q limit, and is converted to P-Q bus, then a new 
row and a new column corresponding to that bus is added to the 


tB’ ’ ] matrix 

in frDL* 

and 

vice versa Since [B’*j 

matrix 

has 

changed , 

it 

has to 

be 

retr langu 1 ated/re inver ted 

this 

will 

increase 

the 

load 

f low 

solution time drastically 

and 

thus 


deteriorates the attractive feature of tBLl* 

in view of the above, a sensitivity method is used m which 
the sensitivity of the P-V buses us determined During the load 
flow solution if the calculated reactive power of a P-V bus 
violates the specified Q limits, then the incremental voltage 
magnitude is determined as a product of the reactive power 
mismatch and the sensitivity to update its specified value instead 
of converting it to P-Q bus Ihe change m specified voltage is 

calculated for a P-V bus 1 as 

AQi 

= S 

v V 

1 


AV 
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where is sensitivity term derived from diagonal element of 
inverse of another matrix identical to tB* 1 but assembled with 
l'-V buses also included Ihis sensitivity matrix is formed only 
once at the beginning of load flow solution and utilised for 
Q-limit check in each iteration Ihus the sensitivity method of 
load flow is much faster 

2 7 SYSlhM SlUDihS 

Various fast decoupled load flow models considering both 
Stott's and Van Amerongen’s assumptions as discussed in previous 
sections CfrJDLt’-l to tDLl*-Vl) have been tested by solving the load 
flows for the following power system examples on H.1' 9000 computer 
system 

1 ihhh 14 bus system 

2 JLfcJbh 30 bus system 

lhe system data is given in Appendix 1 & 2 A convergence 
criterion of 0 0001 p u has been taken using a power base of 100 
MVA lhe results have been obtained for both well behaved and 
l 1 1-conditioned cases ill conditioning in the system have been 
created by increasing 11/X ratio of all the lines simultaneously 
which has been achieved by raising the series resistance of all 
the lines by a common factor ’a' simultaneously from their base 
values 

llesults of unadjusted load flow solutions for the two test 
systems are given in lables 1 & 2 which compares the performance 
of differend load flow models in terms of no of iterations taken 
by each method for both well behaved and ill conditioned system 

lhe sensitivity method have been used to handle Q-limit m 
all the above models and the results of adjusted load flow 



2 * 

solutions comparing no of iterations tak.cn b> the six. versions of 

1 131 t* for only ILfcL If bus sjsltin is given in lable 3 

Approximate CPU time taken b* ea^h iteration of l DLl in 

polar and Jt*DLl in rectangular is 0 08 and 0 0b sec 

respect i veiy 

2 8 Conclusions 

irom the results given in tables (.2 1 to 2 3) we can conclude 

that 

1 lor well behaved system, the number of iteration taken by the 
standard IDLl taking Stott and Alsac’s assumptions is the 
least for base case for both the 14 and 30 bus unadjusted and 
14 bus adjusted cases 

2 lor ill conditioned cases, the General purpose 1-DL1 algorithm 
by Amerongen in polar coordinates takes the least iterations 
followed by the tUbk. V (in hybrid coordinates) 

3 lhe M)SL IV using Amerongen s assumptions in rectangular 
coordinates is the worst for both the 14 and 30 bus systems 
Irom this it is concluded that IDLlr 11 using Amerongen* s 

assumptions in polar coordinates is the best method 
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Comparison of various tDLt 1 models for unadjusted load flow 
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lABLh l 1 


Comparison of various MiLi* models for unadjusted load flow 


solution (JLfc.fc.fc. J0 bus) 
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lABLfc l 3 


Comparison of various tDO* models for adjusted load flow 


solution Htht. 14 bus) 
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NEW tASl DkCOUDLLD SlAlh LSI IMA 1 OH MODLLS 

3 1 INI ItODUL 1 ION 

the problem of monitoring the power flows and voltages on a 
transmission system is very important for maintaining the system 

security By checking each measured value against its limit, 
one can identify where problems exist in the network and can 
take corrective actions to relieve overloaded lines or 

out-of-limit voltages 

Ihese, and many other advanced functions such as 

optimal power despatch, unit commitment, automatic generation 
control etc are performed at the energy control centers and 
requires a reliable system data base Ihe system data acquired 
from field, ma> be corrupted when received at energy control 
centers due to inherent inaccuracies of measuring devices, 
transducers, A/D converters and the telemetry link If the 
measurement errors are small, they may go undetected and can cause 
misinterpretation 

A state estimator can smooth out small random errors in 
meter readings, detect and identify gross measurement errors, 
and fill in the missing readings that might have not been 
acquired due to communication failure One of the state estimation 
methods most popularly being used by utilities at energy 
control centers, is the Jbast Decoupled state estimator ODSh) 

[ ] based on the Weighted Least Square (WLS) technique ihe 
tDSL is quite fast as it uses constant information matrices in 


each iteration However, even 


LDSL method, 


state estimator 
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program may take few minutes to estimate the state of a 
moderate size system whereas the data scan time varies between 

one to ten seconds the research has been proceeding towards 

developing newer and faster algorithms m order to achieve this goal 
but still remains a challenge to researchers in this area 

Narasimhan Vempati et al [34] suggested a new algorithm 
using orthogonal izat ion of the normal equations method to solve 
the static state estimation problem Iheir algorithm is more 
robust than the normal equations method and is claimed to be 
extremely fast m updating the state for small change m the 
measurement set However, they have applied Newton ltaphson (Nil) 
equations for the state estimation Ihe fast decoupled model is 
known to be superior in speed and time over the NK method for 
solving the load flow problem Ihe motivation behind the present 
works was to develop a model based on orthogonal izat ion of fast 
decoupled equations 

A comparison of the various load flow methods shown in the 
previous chapter indicates that the General purpose load flow 
algorithm suggested by Van Amerongen [18] performs better than 
the Stott’s MlLi* [b] specially in solving ill-conditioned cases 
Ihe motivation was to develop IDSJb model using Amerongen’ s 
assumptions and considering normal as well as orthogonal lsed 
equat ions 

Ihus, in this chapter four different models of tDSh has been 
developed, two based on normal equations and the two using 
orthogonal izat ion, utilizing Stott’s as well as Amerongen’ s load 


flow mode l s 
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3 2 Wfc.iGHlfc.D LLAS1 SQUAUfc. Mfc l HOD IN POWfc.lt GY£>lfc.M SlAlfc. EbllHAliON , 

Ihe three criteria on which the static state estimation are 
generally developed are the maximum likelihood criterion, the 
weighted least squares criterion and the minimum variance 
criterion Since the weighted least squares (WLS) criterion is ol 

relevance to this work it is described briefly m this section 
the objective used m the WLS method is to minimize the sum of 
the squares of the weighted deviations of the estimated 

measurements from the actual measurements 

the system data received at the energy control center may be 

i 

i 

corrupted due to various reasons as stated earlier the 1 

difference between the true and the measured quantity can be 
termed as the measurement error and its distribution is 
random We can assume that the probability density distribution 
function of such a random error is a normal (Gaussian) 

distribution 

Let Z mea be the measured value of the system data of 
t ru6 

measurements, Z be the true value and y\ be the random 

measurement of the error Ihe measurement value can be expressed 
as 


z mea 


z true 


+ 


v 


(2 1 ) 


if cr is the standard deviation and a is the variance of the 
measurement error, the objective function can be written as 

2 


min JQO 


N 

m 

--5 

1=1 


[ z mea - z true ) 


(2 2 ) 


Ihe measured quantity could be a function of the state 


variable where h is the functional relationship Ihe measurements 



31 


can be modelled m an unified manner as 

L = h(X) + n (2 J) 

where L is a (mxl) measurement vector, X is (nXl) system state, 
h(X) is the (mxl) vector of non- linear measurement functions and 
Tj is a (mxl) vector of measurement errors 

Assuming sufficient redundancy m the measurement l e ,m>n, 
the best estimate of the state is obtained by minimizing the 
weighted least squares objective function 

mm J(X) = [ Z - h(X) ] if 1 T Z - h(X) 1 (2 4) 


Uewritmg the objective function m an incremental form i 
to accommodate the non-linearity of the problem, the objective can' 
be written as 


mm J(AX) 


A Z - H( AX) 


]* -- 1 [ 


AZ “ E(AX) 


(2 5) 


where Z is the (mxl) measurement mismatch vector U is a 
(mxn)measurement Jacobian X is the (nxl) correction vector to 
the system state and 11 is the (mxm)diagonal covariance matrix 
Setting the gradient of the objective function to zero for mini 
mization, the normal equations method obtains the estimate by 
iteratively solving 


[h 1 U 1 hj AX = H l It 1 AZ (2 b) 

where H.* It * H 1 is known as information matrix 
3 J tASl DhCOUPLliD SI Alt hSUMAlOlt 


in fast decoupled state estimator the measurement vector Z 
is partition into active and reactive parts l e 


z 


f (V,<S) 

p 

2 ; 




g(V,<5) 


Z 


h (V,£) 


( 2 7) 



3J 


where L and L are the sets of active and reactive power 
injections power flow measurement, respectively 

With the above definitions we obtain the Jacobian 


Ht£,V) = 


0 V 

av 





33 

av 



H 11 

H 12 

3 Q 

a q 



H 21 

*22 
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3 \J 

_ 


m 

_ 


2 8 ) 


t*or a practical system, the elements of H^ and H^ are 
generally much larger than those of and H^ respectively 

Hence the latter Jacobian submatrices are neglected 

thus, the eqn (.2 8) takes the form 


H(<5 , V ) = 


U 


11 


0 


0 

H 


22 


(.2 9 ) 


lhe information matrix will be 


1 = 


M 11 U p ti u 


0 


0 


U H 

22 q 22 


12 10 ) 


where it and It are the real and reactive power measurement 
P Q 

covariance matrices respectively Denoting H tl by H p and by 

H , we get two decoupled sets of equations for state estimator 
22 

as follows 


H it H 
p p p 


AS 


a* ii lt At* 


(2 11 ) 



3 J 


u„ it u 
q q q 


6\ 


H" U AQ 

q q 


(2 12 ) 


the active and reactive power injections, tine flows can be 
expressed 

r n n 

p i = v i 2 v j °ij COb c<5 i ~ + 5 v j Bi j blN c<5 ’ _ c * l3) 


, - v] 

1=1 -I 

- “ n 

i = v ! } v j Sj biN l4 . - V - } v j B i j cos - v 

L i=l i=l J 


Q, = 


i* = 

t- j 


Q = 
t J 


(2 14 ) 


V - V V COS (6 - 6 ) 
v v ■) v j 


g - v v b siNta - a ) 

1 J t j V J l J 


(2 15) 


V - V V COS 16 - 6 ) 
t i J i J 


b - V V g S 1 N (6 - 6 ) - y . (2 lb) 

v J l ) \ J i j v ahi j 


where C +jB, is the ljth element of the f v and g +j b is 

i J v J bus i j v j 

the series admittance of the line connected between bus l and j 

lhe elements of the Jacobian corresponding to real power 
measurement (H. ) using fast decoupled assumptions lb] are 


H 

p 


di\ ' 

06~ 

31’v j 

*T - 


(2 17 ) 


lhe elements of [B’l and[B’ ] will be different if the above 
tOSh model is developed based on Amerongen’s assumptions as 

compared to the one developed based on Stott’s assumptions Ihese 
two models have been referred as tDSh-2 and t-DSh-l respectively 
In IDSh-1 model while forming IB’] the series resistances of 
lines, shunts and transformer m phase tappings are omitted, 

while forming [B* ] phase shifters are omitted In IDSJb. -2 model 

while forming [B’l the series resistances of lines are not 

omitted, while forming [B ] the series resistances of the 
lines and phase shifters are omitted 
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Applying these assumptions the fast decoupled state estimator 
model is based to solved to following equations in each 
iteration 


A 6 


h 1 it" 1 H 
r p p 


H 1 U' 1 AZ 

p P 


(2 18) 


AV 


H r It" 1 

<1 


u 

4 


11" 4 AZ 
4 <? 


(2 19) 


3 4 SIA11C SlAlh JhSlJLMAlOlt BAbhD ON OltlHOGONAL lltANStOltMAl ION 

lhe linearized objective function eqn (2 5) can be rewritten 
in terms of square root of inverse of covariance matrix til ) as 


J (AX) = 


i / i 

11- /2 az -it" /2 h AX 


1 1 
U" x2 AZ - U" /2 il AX 


(2 20 ) 


Define A to be an orthogonal tmxm) matrix such that A T A = 1 

using orthogonal transformation based on Given’s rotation 
(Appendix 31), 

T .. 1 

(2 21 ) 


A U" " 2 K = 


0 


where D is a diagonal tnxn) matrix and U is a unit upper 


objective function can be rewritten as 

i 


J (AX) = 


A It" ^AZ - A It" /2 a AX 


orthogonal it} 

of 

A, 

the 

l l 

■ 



AZ - A It H 

AX 

(2 

22) 

Appendix 3 1 

) 

can 

be 


applied to the measurement vector also, 


i 

A It" 2 AZ = 


D" 2 AZ. 


d" 2 az„ 


(2 23) 


Minimization of the objective function leads to 
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U AX = AZ x (2 24) 

Back substitution of L on U provides X and hence the new state 
vector (X) If the Jacobian is not considered to be constant then 
orthogonal izat ion is of both Jacobians and measurement mismatch 
is required in each iteration till convergence is achieved in 

order to save Ll’U time Narasimhan Vempati et el 134] have re line 
arized Jacobian after ever} four iteration 

lhe equation (.3 24) can also be directly used to find the 
change m system state for a small incremental change m 
measurement set 

3 5 1‘UOl‘ObhD Nil W IASI DkCOUBLliD SlAlfc. hbliMAlOll 

In the proposed method, the Jacobian submatrices have been 
made constant by applying fast decoupled assumptions and thus the 
computational time required for repeatedly orthogonal iz mg the 
Jacobian in each iteration can be saved lhe constant real and 
reactive Jacobians are transformed only once to form the 
orthogonal matrices A and A , which are used m successive 
iterations Ihis model has been developed using Stott’s as well 
as Amerongen’s assumptions for [B ] and [B* ’ ] matrices and have 
been referred as iDbh-3 and i Dish- 4 mod el s respectively 

lhe major computational steps for the proposed new fast 
decoupled state estimators are listed below 

1 Head system data including the telemetered measurements 

2 frorm Y l matrix, the constant jacobian submatrices and the 
covariance matrices 

3 Set iteration count X = 0 

4 llotate rows of H into U one by one and form the 

p p 

corresponding orthogonal matrix A 

5 llotate rows of H into U one by one form the corresponding 
orthogonal matrix A^ 



3 / 


b Compute Z to obtain uL and uZ b* applying the 

1 Ip 2 p 

rotation operations of step 4 


7 

I’erform back substitution on A Z on 

ip 

U 

P 

to 

obtain 

AC 

8 

Update bus voltage angles & by A6 





9 

Compute Z^ transform Z to obtain 


and 

Z by 

2q 

applying 


the rotation operations of step 9 





10 

Advance iteration count by K. = 

K + 

1 



11 

Perform back substitution on Z on 

lq 


to 

obtain 

AV 


12 Update V by AV 

13 If convergence is not achieved, go to step - b 


3 b bYSlEM SlUDlfcS 

Ihe four state estimator models tDS»E-l, EDSE-il, 

M)Sfc.-iJLl and tDSh-iV discussed m this chapter have been tested 
on the Ihht 14 and 30 Bus test system on HI’ 9000 computer system A 
convergence criteria of 0 0001 p u on state variables 

(.voltages) has been taken using a power base of 100 MV A A base 
case load flow solution was obtained as the measurement set 

of data required for testing the state estimation programs 

Errors were introduced m both the test systems m the real and 
reactive over injections bus no 4 A b as well as in the real and 
reactive power flows of lines 8 & 15 An error of 3% 

(increase) was considered from the base case values in the real 
power quantities and 2% (decrease) was considered in the reactive 
power quantities Uest of the measurements set were not altered 
for running the state estimation program lhe weight matrix 
elements for the perfect measurements were considered as 1000, 
whereas for the imperfect measurements were taken as 100 lhe 
result obtained from the four methods are presented in table 3 1 
and 3 2 for 14 bus system and table 3 3 and 3 4 for 30 bus 



36 


system iabte 3 1 and 3 3 compares the number of iterations and 
tables 3 A and 3 4 compare the state variables incremental state 
estimation was carried out using all the four models by 
introducing a 5% increase in base case line flow for real power 
line No 10 Using the previously calculated base case state 

variables (Voltage magnitudes and angles), the state estimation 
program was run again to determine the incremental change in 
state variables lhe no of iterations required for calculating 
the incremental change has been given m lable 3 5 

lhe CPU times required per iteration using M)SJb. based on normal 
equations are 0 04sec and 0 047sec for 14 and 30 bus system 

respectively and using orthogonal methods are and sec 

respect ive ly 
3 5 Conclusion 

Prom the results presented in table 3 1 to 3 5, one can 

conclude the following 

l) lor the full state estimation run the number of iterations 
required by frl)Si.~2 and M)Sb-4 (both based on Amerongen’s 
assumptions) take the minimum iterations compared to other two 
methods lhe orthogonal method using Stott's assumptions (tDSk-3) 
takes two iteration less as compared to the normal equations 
(Mibii-l) method for 14 bus system However, they take same number 
of iterations for 30 bus system 

ii ) tor an incremental change m the measurement set, the t-DSh-3 
performs worst m case of 14 bus system, whereas all other methods 
require identical number of iterations for both 14 and 30 bus test 
systems 

m) A comparison of Cl’U time reveals that orthogonal method, as 
formulated m the present work, takes slightly more time per 
iteration compared to the normal equation method, which can be 
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improved by incorporating the enhancements suggested in ref 
and the sparsity techniques 

in view of the above observations it is recommended to 
FUSh-4 model 
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use 














or-^^ncr**u^ocr»cDnoc>coiA 
0(h(UN<0000^^0'<00'tOlft 
o^nmoncor-n^ocuor- ^ 

ovv-iftooflo^t^ooncDn ^ 

oorroDO^O' vj 

owoooo nn^naji-ojTr ~ 

oowwww wf^wwoiw w r^ 


r^\ 

t 


<j) 

& 


i i i ! i i i } i i i r i 
n a u « ii n n if ii (i u a a u 


eeeseeEeeeeeee 

lAWWUltAIDIAIAIAQIlfllflUIUI 

tf) 04 >qitl> 4 > 4 » 4 » 4 >ail»&<tJ 4 » 


XIXl'tJ'U-aXJTJXIXinCJTJTJTJ-D 


<^inm^rt7»r»oo-‘Oii^mr-r^<» 

(Mn<M^Oh-t^O'Q«£>00'vD 

^ftini^inr-oojr-^com^N 

wuu^-^^tpjcunnnn^ 

oooooooooooooo 


t 


J 
< 
2 
\ J 
vr 

o 

X 


0 


Tf; 


« a a ii ii » » a n a a » a a 


ui 

uJ 

IU 

H 

Oi 


* 

K) 

Ui 


J 

CQ 


uj 

w 

J 

cfl 

< 

5 

< 

> 


a) 

h 

< 

h 

V) 


it 

o 


"2 

O 

u\ 

W 

< 

CL 

£ 

<g 


aeeaaeaaeasEFS 

muiwininu)utu)iAttiu>wtnu> 

a*«u<ut^^ 6 a>*cidiCL»ai 4 itsi 


>>>>>>>>>>>>>> 



oM < unff'onftiM < io»-to 

Oh-NOMJ&oWtothCOCO^ 

I I I I 1 t I I t I I I I 

a ii ii tr ii ii a n a « a a n a 


6 E 6 B 66 £ 6 eS 6 eee 

taminuiwiniDuiuiuiiainuiin 

000ty4i&df004f000dJ 


V W W W V« W W v w *— * w w 

QOOOQaOQOOQQOQ 


ononw(\i(u>-^^to^co 
o(oncoa>^nNnwo^t. ui 

oh-oooonn^JMcu-»-cu'T 

OOWIU W«JCU-^ W(U (Uft/tU 

I I I I I I I I I I I I 1 

H H II » 0 » « H H II H II 0 If 


P 

J 


£ 


Z 

a 


< 


<r 

Ui 


<C 

£ 

cz 

o 


i. 


ee eeessesEeEee 

'K' w w u W W W w Vw W W W w w 

OOQQaOQOOQOQOQ 


otrotWMo^o^^Mum 

6 uiw^**ta'OTrtto^oo 

oooooooooooooo 


II » (I 0 II ii I! if 0 (i 0 0 


II 


E££Ee£g££g£££g 

mia«ui<ftw«ui^u»ww0iw 

^ di ftf U) ill Hi (t> a 0 4> 4) OP 0 <1 


>*"• "W- W V* w W -w* >«*■ 'm. w Wr v W 

>>>>>>>>>>-^*>>> 


oocoot^^ a m r t 1 i- 

o u> o o o *u ui r~ - o o * o t*> 

y) t ^ ui m oj ui to ti i to f ( 

OOOOOOOOOOOOOO 


if 0 ii 0 it il 0 i U M I! 11 ti 1 

njt »o uJ <7 ^ m 0 ? 4 » is 4 ij m 8 m rj 

I 4* 0 01 0 «U 0 0 Of Ii ty $ « ti £ 

EEESESEF^prrg-r 

>>>>>>>>>>>>>> 


tr 

*2 

O 

P 

a. 

£ 

kA 

\J) 

<c 

t 

‘J 

p 

V/l 

vjr 

2 

? 

3 


r-tun^invoNco^or-aintr 

» n ii » n n h n u » » b n n 
a)<o<o<ou><oa) 07 oito«otocoto 
r=>»=>3DZ)i=)D33^r)3D 
aocQOQCQOiCQaidQmoQffiaimoQ 


Si X X 
ISISI^IO 
Si Si t> CD 
S^Si> 
Si Si X 
Si X Si SI 
Si Si N N 

i I f 
U II It II 
^ 

+J 4 J +J 4 -) 

a b a a 

WWW w 
0 0 0 0 
M/ w w 

-d -d TJ T 3 


^1 OJ *J< 3) t> ** '“'I X 

CO 20 03 s 03 ^ 9 H ^ H 

33X«t sO ^ ^ 

ISISI^ ^ J3-7NHM0* 
*N3 ^ 'NJ 1-4 r-» 'N 'M "N *M Ni 

III f I I I I I I 

u u u ii ii n ii n n ii 

^ S~\ /^N /"\ <*\ ^*>. <^N A 

4J+J4J V +J P 4J 4J +J 4J 

s s a a a a a a a a 

CO c /3 W CO H) ffl tfi tfl W tQ 

000 000000)0 
\^f W W 

-a*d , d‘d*d*dTJ , d*d*d 


s' 40 


IU 


N Ifl 33 t' f' 5130 

x x x ^ t- t- «* 

t> X X *i* Si N X 

^ h^)33 

•M N ^ X X X -h 

Si Si Si Si Si Si Si 


O) 9 ^ N tf 0 

CO to N ^ 0 SI 0 

X t- t- ^ ^ -ri lO 

3 ^ t- 'M 'O SI 

O) N °0 30 W h 9 

lO to to W c^- 

Si Si SI Si Si Si Si 

,4 H H r4 H t4 tH H H H H rl H 

ii ii it li ii ll ii H ii I' I' 

^ s~\ s~\ s~\ ^ ^ 

4J +J 4J +J +J 4J 4J 

sssBs^ssaasaae 
MtorawwMM mtowMwrow 

0000000 0O00000 

^ ^ ^ *^/ \^/ v v</ •«»/ v v v-> v 

>>>>>>> >>>>>>> 


Si-^oxSixm'M'MSic-of-^ 

SiXX^XCO»-it^CDC^'-iX' r OX 

Sit>c-Sit-t>^iH'V)SiO)t-xSi 

Sj^hMhhNhh^wthhM 

i 1 I l i I f I i i i f f 

U li H H (I ft I! II || II H fl tl 11 

SSSeSSSEESSSSES 

Pf 3 +J 4 J^J 4 J 4 J 4 J+J+J 4 JPfl-P 

W ffl M tfl B tn B } to to to CO (0 CO D 9 

00000000000000 
%w/ V w W W W S-/ W Sw/ W w *S^ w 

OQQQQaqooqqqqQ 


SlXSiX^'MN^i-433'tfSiO)X 

SixxxcDdO'Oorj'NJSix^m 

§£*SiSlSiSi'^' , OX'^N’H'M^ 

SUSNNNMNhh'NNM^^ 

I ! 1 I t I i I i I ( I I 

U It H I! If II II II |( II || || || || 

^ /“N /~\ /*\ 

cdttJdoJcJfljflJctJcdctJaJaJedd 

00000000000000 

aaaEESEsaaaEEE 

w/w^s-,/WV^W 


SitrtriilOijQXCO'MCOXXSlXSi 

Qin^^'Mioco'^mo'NtjosiJD 

x^sO^xiOTfn'M^iotnin^ 

SlSiSiSiSiSiSiSlSiSiQSiSlSi 


II II II I! tl II ll II II I! II II tl II 

f~s S*> »"*‘\ S**s ^ /*> ^ /--N 

aEassssasaesEB 

+JP+J+J 4 J 4 J 4 JP 4 J+ 5 + 04 J+J 4 J 
0100000030000000 
00000000000000 
v-/ v/ <«/ w w W 

>>>>>>>>>>>>>> 


S!SlSiSl 5 J^ , '<*XvQ' v )‘r)rfi^*H 

SitnSiSiSi^iot^^SiviD^Sitn 

0 tl»CsHO 300 NN 0 tO 0 im 

SiSiSiSiSiSiSiSlSiSiSiSiSlSi 


it li it ti ll it il li li ll u tl ii u 

'•“>. •'*> .**> /> ^"S /**\ <<■>. v^> /*\ S*\ /N 

ddtdotfdccJrfocJBjccJccJtfStdKJ 

00000000000000 

aaaeaaaaaaaaaa 

V W/ VW w/ vw w W W V 
>>>>>>>>>>>>>> 


HN'7^in0C-0O39H'M'O4j) 

II II II >1 11 II U li ll II ll 11 II II 
A'AJiJi'AJiJijyjyjiJi'A&ji 
2DDDDDDDDDDD32 
00000000000000 


£ 


j 

< 

o 

O' 

O 

X 

b 

o 

v7 


0 

e 

1 \ 

X 


2 

b 

J> 

‘<ar 

cu 

J 

< 

E 


H 

2. 

Cj 

/- 

cy, 

E 

VI 

<c 


2 

JA 

V5“ 

2 

o 

ul 

y 


vr* 

z 

V/i 















T/mu: i 4 


4 1 


Co^PA^I^CjN S-TAfE \MR.i AR>l(~ s ftt t tvitc i-J w » T » s, AMD 

A^E£.Ov4^r.wi ^ErHOt> BY CftTHo^OHAL TECHtfjqc. 

(_s.0 RviX 'TEi.r system} 


ST^TTs* method 


AvAE-dOK^t w!s Piernoii 


BUS 8 ? I 
BUS* £ 
BUS® 3 
BU=>* 4 
BUS® 5 
BUS® 6 
BUS® 7 
BUS- 3 
BUS* 9 
BUS® 1 0 
BUS® ? 1 
8U c *te 
BUS® t 3 
BUS=? *1 
BUS® \ $ 
BUSat 6 
BUS®! 7 
BUS®! 3 
BUS® J 9 
BUS®C 0 
BUS =2 1 
BUS=BS 
BUS«£3 
BUS®£4 
BUS®25 
BU5«£6 
BUS *2? 
BUS®23 
BUS*£9 
BUS*30 


VCestm.)®! 0600 DCestmJ® 0000 
V I est m ) » 1 0457 D(estm) = - 0752 
V{ estm)»1 0159 D(*stm)=- 1404 
VCestm)®! 06401 DC estm ) = - 2021 
V< ! 0 J 2 1 Dtestm )*- 2156 
VCastm)*! 0537 D<*stm)*- 2202 
VI<»8tm)*1 0054 D( ©atm >**- 1842 
V<e*tm)«l 0245 D(estro)*- 1063 
VCeatm)»1 0417 B(estm)*- £011 
VI es t m ) * 1 0366 Dlestm)*- 2288 
VIestm}*! 0158 Dtestm)*- 1275 
VI**tm>®1 0474 DC estm }*~ 2201 
Vteetm)*! 0145 DCestmJ*- 1458 
VIestro)*! 031? Dtestm)*- 2341 
VI*stm>*1 0273 2343 
V ( est ro )*1 0345 DIestm) 88 - 2285 
Vleetm)*! 030? D(estm)*- 2325 
VIe«.tm)«l 0156 DIestm)®- 2442 
V(estm)®} 0132 DIestm)*- 2470 
Vlestro)*! 0186 D I est m ) =~ 2438 
VIestm)*» 0237 DIestm)®- 2351 
VI estm)=1 0244 D( estm)=- 2344 
VI estm)*1 0193 DIestm)*- 2358 
VIestm)*! 0188 DIestm)*- 2317 
VI estm)»l 0287 DIestm)*- 2016 
V(estm)»l 022? DIestm)*- 2012 
Vleatm)® 1 0352 DIestm)*- 1874 
VIestm)*! 0141 D(estm)*- 1490 
VIestm)*! 0261 D( ©stm)=- 1925 
VIestm)*! 0181 DIestm)*- 2028 


VCesr»«i)*l 0600 D(estra)® 0000 
VI est m ) * 1 0456 DCestm )*- 0751 
V ( sstm ) *1 0156 DIestm)®- 1404 
Vieatm)®1 0639 DCestm)®- 2025 

V c e stm ) *1 0120 DCestm)®- 2155 
Vreatm) - ! 0586 D ( est m ) =~ 2195 

0 OSS DCestm)®- 1843 
Vf estrn ) ~1 0243 DCestm)®- 1064 
VCestm)*1 0414 DCestm)®- 2008 
V(estra )=1 0363 DCestm)®- 2285 
vr estm )* 1 0155 DCestm)®- 1278 
VCestm)®? 0473 DCestm)®- 2191 
Vi*3tm?®1 0143 DCestm)®- 1460 
VCe$tm)*1 0307 D( estm)*- 2333 
VCestm»-1 026"’ DCestm)*- 2336 
VCestro)*1 0341 DC ?stm )«- 2281 
Vtestm)®! 0297 DCestm)®- 2322 
VI© 3t«)*1 0152 DCestm)®- 2439 
Vf&stm>*1 0127 DCestm)®- 2469 
Vl©5tm)*1 0180 DCestm)®- 2436 
Vfssn?i)*l 0238 DCestm)®- 2350 
Vfe«sf‘it >*1 023« DCestm)®- 2342 

V < e s t m ) - 1 0188 DCestm)®- 2351 
VIestm)*! 0183 DC estro )=- 2306 
Vf*»9tm}®1 0282 DCestm)®- 1995 
Vf estm ) *1 0218 DCestm)®- 1981 

V f *stm ) =1 0349 DCestm)*- 1854 
Vfesfm)*i 0140 DCestm)®- 1491 
VCestm)®! 0258 DCestm)®- 1898 
vc estm )*1 0181 DCestm)®- 1993 
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Comparison of various IDSL methods for incremental change 


IDSL Methods 

No of iterations 

14 bus 

30 bus 

IDSL 1 

4 

z 

IDSL II 

2 

A 

IDSL III 

4 

2 

IDSL IV 

2 

2 
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CONCLUSION 

lhe work carried out in this thesis was aimed at comparing 
different fast decoupled load flow methods in polar and 
rectangular coordinates as well as exploring new solution 
techniques for power system load flow studies lhe standard JbDLt 
by Stott and Alsac does not behave well in ill conditioned cases, 
whereas the General purpose IDLI by Amerongen can handle even if 
the network system is ill conditioned lhe potential of existing 
IDLI methods using Stott’s as well as Amerongen s assumptions have 
been compared with the new algorithms on Ifc.fc.fc. 14 and 30 bus 
systems for both well and ill conditioned cases 

A comparative study of the performance of the various fcULfc 
techniques indicated the best method which could be used to 
develop a new fast decoupled state estimator 

Narasihman Vempati et al [341 have applied the 
orthogonal isat ion technique to the state estimator However their 
work has been in Newton method Ihey have claimed that the 
orthogonal isat ion makes the algorithm more robust in ill 
conditioned cases in this thesis the orthogonal isat ion by Givens 
rotation has been applied to the fast decoupled state estimator 
using Stott's as well as Amerongen’s assumptions Hence, an 
attempt has been made to explore the possibility of combining 


both the advantages 



frrom the comparison of the test results of the various IDLl 
and 1 DSh methods presented in lables 2 1 to J 5 , it is concluded 
that 

1 lhe general purpose fast decoupled load flow by Van Amerongen 
is found to take the least iterations m ill conditioned systems, 
however the standard IDLl* by Stott and Alsac takes the least 
number of iterations for well behaved systems 

2 Between the two algorithms using Amerongen’ s assumptions, the 
one m polar coordinates CtULt II) takes lesser number of 
iterations than IDLI VI 

J lhe time taken per iteration by the IDL! methods m polar 
coordinates is slightly greater than m rectangular coordinates 

4 lor the base case study of fast decoupled state estimator, 
the method using orthogonal l sat ion by Givens rotation taking 
Amerongen s assumptions is the best for both the base case and 
incremental estimation lhe ID&h method by normal equations using 
Stott's assumptions is the worst in lfc.fc.li 14 bus test systems 
however all the methods perform equally well for the 30 bus test 
system 

5 lhe time taken by the state estimator using normal equations 
is lesser than using the orthogonal isat ion method 

lhe orthogonal isat ion procedure is computationally intensive 
and hence the time taken by this method is more than the normal 
equations method Vempati et al [34] have however described the 
enhancement to the Givens rotation by which the number of 



4b 


multiplications is reduced they have also exploited the sparsity 
of the Jacobian by applying various sparsity techniques which has 
not been done in this thesis 

lhe bad data measurement is one of the foremost problems to 
be tackled m the study of state estimation Narasihman Vempati 
et al [34 1 have claimed that by applying the orthogonal isat ion 
principle to the state estimator problem, the identification of 
bad data and its removal is easier Hence, the present work may be 
extended in future to incorporate this important feature in fast 
decoupled state estimator programs, which could enhance its 
utility 
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A1T1ND1X ft 

Consider the following two row vectors 

U = [ 0 , , , 0 ,U v , , , U l , , , U f +i j 
K = [ B , . , B , X,, , j 

A plane rotation boLwcrn tl and X is dc 1 me d sui li Lhat the 
ith element of X is annihilated After the rotation the row vectors 
take the form 


u-= [ 0 , , , a ,u;, , , u;, , , ] 

*’= [ 0 • - ■ 0 ' a - • ■ V ■ ■ ] 


the two vectors U and X are defined as 



C is ' 


' u ' 


u’ ■ 


-is C 


X 


X’ 

where 

2 2 

C + is = 

1 





Ihe scalars C and is are 
determined from the requirement that X^= 0 and are 
given b$ 

X. 

C = 

^ u 2 + x 2 

X 

1. 

is = — 

** u 2 + x 2 

•L C 

lhe above is known as orthogonal transformation using Civens 


rotation 
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DATA FOR IEEE-14- BUS TEST SYSTEM CAT 100 MVA BASE} 


j3-1 Bus Data 


Bus No 

Generation 

Load 


HR FTOAR 

MW MVAR 


1 

0 0 

0 0 

0.0 

0 0 

S 

40.0 

0 0 

21 7 

12 7 

3 

20 0 

0 0 

94 H 

19 0 

4 

0 0 

0 0 

11.2 

7 5 

5 

0.0 

0 0 

0.0 

0 0 

6 

0 0 

0 0 

47 8 

-3 9 

7 

0.0 

0 0 

0.0 

0 0 

8 

0.0 

0 0 

7.6 

1 6 

9 

0 0 

0 0 

29.5 

16 6 

10 

0.0 

0 0 

9 0 

5 8 

11 

0.0 

0.0 

3.5 

1 8 

12 

0,0 

0.0 
'o 0 

6.1 

1,6 

13 

0.0 

1J.5 

5.H 

14 

0.0 

0 0 

14.9 

5 0 



B-F Line Data 
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Line No 

Pub on 
From 

To 

Line Impedance 

R in p u X in p u 

Half Line 
Charging 

sus ceptance 
( p. u) 

1 

1 

2 

0 01938 

0 05917 

0 0264 

E 

1 

8 

0 03403 

0.22304 

0 0246 

3 

2 

3 

0 04699 

0 19797 

0 0219 

4 

2 

6 

0 038 11 

0. 17632 

0 0187 

5 

2 

8 

0 05695 

0 17388 

0 0700 

6 

3 

6 

0 06701 

0 17103 

0 0173 

7 

4 

11 

0 09498 

0 19890 

0 0 

a 

4 

12 

0 12291 

0 25581 

0 0 

9 

4 

13 

0 06615 

0 13027 

0 0 

’10 

6 

7 

0 0 

0 20912 

0 0 


6 

8 

0 01335 

0.0421 1 

0 0064 

12 

6 

9 

0 0 

0.55618 

0 0 

13 

7 

5 

0 0 

0.17615 

0 0 

14 

7 

9 

0 0 

0.11001 

0 0 

13 

8 

4 

0 0 

0.25202 

0 0 

16 

9 

10 

0 03181 

0 084^0 

0 0 

17 

9 

14 

0 12711 

0 27038 

0 0 

18 

10 

11 

0 08205 

0 19207 

0 0 

19 

12 

13 

0 22092 

0.19938 

0.0 

20 

13 

14 

' 0 17093 

0.34802 

0.0 





B-3 

Transformer Data 


Tram f ormer 


Buses 

Tap setting 

From 

To 

1 

6 

7 

0 978 

2 

6 

9 

0 969 

3 

8 

4 

0.932 


ft-4 Shunt Capacitor Data 


Bus No. 

Susceptance p u 

9 

0.190 



B-5 Voltage-Controlled Bun Data (P-V buses) 


Bus No. 

Vo 1 tage 
magnitude p u 

Reactive 
Minimum MVAH 

Power Limits 

Ma<imuni RVaR 

2 

1.045 

- 40 0 

50 0 

3 

1 010 

0 0 

40.0 

4 

1.070 

- A 0 

?4 0 

5 

090 

- 6 0 

24 0 


Slack bus voltage - 1 06/0° 
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APPENDIX - C 


DATA FOR IECE-30 BUS TE9T SYSTEM (AT 100 MV A BASE} 


C- 1 Bus Data 


Bus No. 


Generation Load 

HE FTOAR MW MVAR 


1 

0 0 

0 

0 

0 

0 

0 . 

0 

2 

40.0 

0 

0 

21 

7 

12 

7 

3 

0 0 

0 

0 

30 

0 

30 

0 

4 

0 0 

0 

0 

0 

0 

0 

0 

5 

0 0 

0 

0 

94 

2 

19 

0 

6 

0.0 

0 

0 

0 

0 

0 

0 

7 

0 0 

0 

0 

E 2 

a 

10 

9 

a 

0 0 

0 

0 

2 

4 

1 

2 

9 

0 0 

0 

0 

0 

0 

0 

0 

10 

0 0 

0 

0 

5 

8 

2 

0 

ii 

0 0 

0 

0 

7 

6 

1 

6 

IE 

0 0 

0 

0 

11 

2 

7 

5 

13 

0 0 

0 , 

.0 

0 

0 

0 

0 

14 

0 0 

0 

0 

6 

2 

1 . 

.6 

15 

0.0 

0 

0 

a 

n 

C. 

2 

5 

16 

0.0 

0 

0 

3 

5 

1 

8 

17 

0 0 

0 

0 

9. 

.0 

5 

.8 

IS 

0 0 

0 

0 

3 

2 

0 , 

.9 

19 

0.0 

0 

0 

9 

5 

3 , 

.4 

SO 

0.0 

0 

0 

2 

2 

0 . 

.7 

21 

0.0 

0 

0 

17 

5 

1 1 . 

.2 

EE 

0.0 

0 

0 

0 

0 

0 

0 

23 

0.0 

0 

0 

3 

2 . 

1 , 

.6 

24 

0.0 

0 

0 

a , 

.7 

6 

7 

25 

0.0 

0 

0 

0 , 

.0 

0 , 

.0 

26 

0 0 

0 

0 

3 

5 

2 

3 

27 

0.0 

0 

0 

0 

0 

0 

0 

28 

0.0 

0 

0 

0 

0 

0 , 

.0 

29 

0.0 

0 

0 

2 . 

.4 

0 , 

.9 

30 

0.0 

0 

0 

10 

6 

1 , 

.9 
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C-p Line Data 


Dugas Lina Impedanrp HMf Line 

Line No. From To R in p.u X 'in p.u. Charging 

sus ceptance 
( p u) 


1 

1 

2 

0 0192 

0.0575 

0 0264 

2 

1 

8 

0 0452 

0. 1852 

0 0204 

3 

2 

11 

0 0570 

0. 1737 

0 0 184 

4 

8 

11 

0.0132 

0.0379 

0 004 

3 

2 

5 

0 0472 

0 1983 

0 020 y 

6 

2 

13 

0 0581 

0 1763 

0 0187 

7> 

11 

13 

0 0119 

0 0414 

0 0045 

8 

5 

7 

0 0460 

0 1160 

0 0102 

9 

13 

7 

0 0267 

0 0820 

0 0085 

10 

13 

3 

0 0120 

0.0420 

0 0045 

11 

13 

9 

0 0 

0.2080 

0 0 

12 

13 

10 

0 0 

0 5560 

0 0 

13 

9 

4 

0 0 

0 2080 

0 0 

14 

9 

10 

0 0 

0 1100 

0.0 

15 

11 

12 

0 0 

0 2560 

0 0 

16 

12 

6 

0.0 

0 1400 

0 0 

17 

12 

14 

0.1231 

0 2559 

0 0 

18 

12 

15 

,0 066? 

0 1304 

0 0 

19 

12 

16 

0.0945 

0 1987 

0 0 

20 

14 

15 

0 2210 

0. 1997 

0 0 

21 

16 

17 

0 0824 

0. W 23 

0 0 

22 

13 

18 

0 1070 

0 2185 

0 0 

23 

18 

• 19 

0 0639 

0.1292 

0 0 

24 

19 

20 

0 0340 

0.0680 

0 0 

23 

10 

20 

0 0936 

0.2090 

0 0 

26 

10 

17 

0 0324 

0 0845 

0.0 

27 

10 

21 

0 0348 

0 0749 

0 0 

28 

10 

22 

0 0727 

0 1499 

0.0 

29 

21 

22 

0 0116 

0 0236 

0.0 

30 

15 

23 

0.1000 

0 2020 

0 0 

31 

22 

24 

0 1150 

0.1790 

0.0 

32 

23 

24 

0. 1320 

0 2700 

0 0 

33 

24 

25 

0. 1885 

0.3292 

0 0 

34 

25 

26 

0 2544 

0.3800 

0 0 

33 

23 

27 

0 1093 

0.2087 

0 0 

36 

28 

27 

0.0 

0.3960 

0.0 

37 

27 

29 

0.2193 

0 4133 

0.0 

38 

27 

30 

0 3202 

0.6027 

0 0 

39 

29 

30 

0 2399 

0 4533 

0 0 

40 

3 

28 

0 0636 

0.2000 

0 0214 

41 

13 

28 

0.0169 

0.0599 

0.0063 
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C-3 Transformer L'nta 


Transformer 


Luses 

Tap setting 

From 

TO 

1 

11 

12 

0.932 

2 

13 

9 

0 978 

3 

13 

10 

0 969 

4 

28 

27 

0.968 


C-4 Shunt Capacitor Data 


Bus No. Susceptance p.u 


10 0 1^0 

24 0 043 


C-5 Voltage-Controlled Bus Data (P~V buses) 


Bus No. 

Voltage 
magnitude p u. 

React t ve 
Minimum MVAR 

Pouer Limits 

Ha < i mum MVAR 

H 

1.045 

- 40 0 

50.0 

3 

1.010 

- 10 0 

o 

o 

4 

1.002 

- 6 0 

P4.0 

5 

1,010 

- 40.0 

40.0 

6 

1.073 

- 6 0 

24 0 


Slack bus voltage = 1.06/0° 
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